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ABSTRACT 

We present four new exoplanets: HIP 14810 b & c, HD 154345 b, and HD 
187123 c. The two planets orbiting HIP 14810, from the N2K project, have 
masses of 3.9 and 0.76 Mj^p. We have searched the radial velocity time series 
of 90 known exoplanet systems and found new residual trends due to additional, 
long period companions. Two stars known to host one exoplanet have sufficient 
curvature in the residuals to a one planet fit to constrain the minimum mass of the 
outer companion to be substellar: HD 68988 c with 8 Mj^p < m sin i < 20 Mj^p 
and HD 187123 c with 3Mj^p < msini < 7Mj^p, both with P > 8 y. We have 
also searched the velocity residuals of known exoplanet systems for prospective 
low-amplitude exoplanets and present some candidates. We discuss techniques for 
constraining the mass and period of exoplanets in such cases, and for quantifying 
the significance of weak RV signals. We also present two substellar companions 
with incomplete orbits and periods longer than 8 y: HD 24040 b and HD 154345 
b with msini < 20Mj^p and msinz < 10 Mj^p, respectively. 
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Introduction 



Of the 151 nearby stars known to harbor one or more planets, 19 are well-characterized 
multiple-planet systems, a nd an additional 2 4 show radial v elocity (RV) residu als indicative 
of additional companions (IButler et al.ll2006l ). For instance, IVogt et al.l (120051 ) reported ad- 
ditional companions around five stars, including two revealed by incomplete orbits apparent 
in the RV residuals (HD 50499 and HD 217107), and one as a short-p eriod, low arn p litude 
variation in the residuals of the fit to a long-period outer companion. iRivera et al.l (120051 ) 
detected a 7.5 Earth-mass mass companion to GJ 876 in a 2-day period through analy- 
sis of the RV residuals to a 2-plan et dynamical fit of the more massive, outer exoplanets. 
Gozdziewski fc Migaszewskil (120061 ) similarly analyzed the RV residuals of 4 stars to search 
for Neptune-mass companions. 

Very little is kn own about the frequ ency or nature of exoplanets with orbital distances 
greater than 5 AU (IMarcy et al.l l2005bl ). Precise radial velociti es have only reac hed the 
precision required to detect such objects within the last 10 years (IButler et al.lll996l ). which 
is less than the orbital period of such objects (P > 12 y for exoplanets orbiting solar mass 
stars). Thus, the RV curves for such planets are all necessarily incomplete, and we must 
obtain many more years of data before our knowledge of their orbits improves significantly. 

The ability to put constraints on planets with incomplete orbits, however weak, allows 
us to peek beyond the 5 AU completeness limit inherent in the ten-year-old planet searches. 
Characterizing incomplete orbits also increases our sample of known multiple exoplanetary 
systems, which improves our understanding of the frequency of orbital resonances, the growth 
of multiple planets, and the mechanics of orbital migration. 



In this work, we present our analysis of the RV data of iButler et al.l (120061 ) in an effort 
to determine which of those systems have additional, low-amplitude companions. 

Many systems known to host one exoplanet show more distant, long-period companions 
with highly significant but incomplete orbits. In these systems, it can be extremely difficult 
to constrain the properties of the outer companion: in the case of a simple trend with no cur- 
vature, very little can be said about the nature of these companions beyond their existence, 
but even this informs studies of exoplanet multiplicity and the frequency of exoplanets in 
binary systems. 

In §[2] we discuss a new multiple planet system from the N2K project, HIP 14810. In §[3] 
we describe how we have employed a false alarm probability statistic to test the significance of 
trends in the RV data of stars already known to host exoplanets. We find that six stars known 
to host exoplanets have previously undetected trends, and thus additional companions. 
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When the RV residuals to a single Keplerian show significant curvature, one may be able 
to place additional constraints on the maximum m sin i of the additional companion. In §§ H]- 
E]we present our analysis of this problem in the cases of HD 24040 b and HD 154345 b, two 
substellar companions new to this work with very incomplete orbits. By mapping space for 
Keplerian fits, we show that HD 154345 b is almost certainly planetary (msini < lOMj^p), 
and that HD 24040 b may be planetary (5Mj^p < msini < 30Mj^p). 

In § we describe how we extended this method to the RV residuals of known planet- 
bearing stars which show trends. We find that for 2 stars we can place sufficiently strong 
upper limits on msini to suggest that the additional companions are planetary in nature. 



2. HIP 14810 



HIP 14810 is a metal-rich ([Fe/H] = 0.23) G5 V, V=8.5 star wh ich we have observed 
at Keck Observatory as part of the N2K program (IFischer et al.ll2005l ) since Nov 2005. Ta- 
ble [1] contains the RV data for this star. Its stellar characteristics are listed in Table [2l 
determined using the sam e LTE spectral analysis used for stars in the SPOCS catalog 
(jValenti fc Fischerl l2005al ). We quickly detected a short-period, high-amplitude compan- 
ion (P = 6.67 d, msini = 3.9Mj^p) and a strong, ~ 200 m/s trend. Further observations 
revealed evidence for substantial curvature in the residuals to a planet plus trend fit. FigH] 
shows the RV curve for this star decomposed into Keplerian curves for the b and c compo- 
nents, and Table [3] contains the best-fit double Keplerian elements. 
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Table 1. RV Data for HIP 14810 



Time 


Radial Velocity 


Unc. 


(^JJJ-z44UUUU ) 


(m/s) 


(m/s) 


13693.760579 


-130.8 


1.3 


13694.831481 


-473.6 


1.2 


13695.909225 


-226.9 


1.2 


13723.786250 


162.6 


1.0 


13724.688484 


324.9 


1.2 


13746.814595 


2.4 


1.3 


13747.852940 


-435.82 


0.94 


13748.734190 


-433.3 


1.2 


13749.739236 


-71.4 


1.2 


13751.898252 


358.3 


1.1 


13752.807431 


241.05 


0.80 


13752.912477 


211.6 


1.7 


13753.691574 


-79.8 


1.1 


13753.810359 


-137.6 


1.2 


13753.901042 


-180.2 


1.2 


13775.836157 


-240.01 


0.97 


13776.812859 


123.4 


1.4 


13777.723102 


346.7 


1.3 


13778.720799 


416.0 


1.3 


13779.744410 


238.4 


1.3 


13841.722049 


-515.7 


1.4 


13961.130301 


-280.9 


1.0 


13962.133333 


-413.4 


1.1 


13969.097315 


-348.3 


1.2 


13981.969815 


-476.5 


1.2 


13982.947431 


-200.9 


1.2 


13983.981470 


151.5 


1.0 


13984.096979 


187.3 


1.2 


13984.985775 


345.7 


1.3 


13985.102106 


357.6 


1.3 
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A 2-planet Keplerian fit yields an outer planet with m sin i = 0.95 Mj^p, P = 114 d, and 
eccentricity of 0.27. We present the orbital solutions for this two-planet system in Table [31 



3. Detecting Long-Period Companions 

Very long period substellar companions appear in radial velocity data first as linear 
trends (constant accelerations), then as trends with curvature, and finally, as the duration of 
the observations becomes a substantial fraction of the orbital period, as recognizable portions 
of a Keplerian velocity curve. It is important, then, to have a statistically robust test for 
trends in velocity residuals. In this section, we discuss calculating false alarm probabilities 
(FAPs) for such trends. 



3.1. Using FAP to Detect Trends 



Marcy et al.l (§ 5.2 l2005bl ) present a detailed discussion of using false alarm probabilities 
(FAPs) for determining the significance of a periodic signal in an RV time series. Here, our 
task is similar. We wish to test the hypothesis that a star has an additional companion with 
a long period, manifest only as a linear trend in the RV series. We compare this hypothesis 
to the null hypothesis that the data are adequately described only by the best-fit Keplerians 
and noise. 

We first fit the data set with a Keplerian model and compare the xl statistic to that 
of a model employing a Keplerian plus a linear trend. If this statistic improves, that is, if 

~ trend -^^no trend iisgative, then the inclusion of the trend may be justified. 
To test the significance of the reduction in x^, we employ an FAP test. 

We first employ a bootstrap method to determine our measurement uncertainties. We 
subtract the best-fit Keplerian RV curve from the data and assume the null hypothesis — 
namely that the residuals to this fit are properly characterized as noise and thus approximate 
the underlying probability distribution function of the noise in the measurements. We then 
draw from this set of residuals (with replacement) a mock set of residuals with the same 
temporal spacing as the original set. 

By adding these mock residuals to the best-fit Keplerian RV curve we produce a mock 
data set with the same temporal sampling as the original data set, but with the velocity 
residuals "scrambled" ("re-drawn" might be a better term since we have drawn residuals 
with replacement.) It is important in this procedure that internal errors remain associated 
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Fig. 1. — RV curve for HIP 14810 with data from Keck, showing the inner planet with 
P — 6.67 d and msini = 3.9Mj^p the outer planet with P — 95.3 d and msini = 0.76 

Mjup- 
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with the scrambled residuals. This ensures that points with error bars so large that they 
contribute little to the xl sum, but nonetheless lie far from the best-fit curve, do not gain 
significance when "scrambled", inappropriately increasing xl- 

We then compare Ax'i for our mock data set to that of our genuine data. By repeating 
this procedure 400 times, we produce 400 mock sets of residuals and 400 values for Axl- If 
the linear trend is simply an artifact of the noise, then re-drawing the residuals should not 
systematically improve or worsen Axl- Conversely, if a linear trend is significant, then the 
null hypothesis, that the residuals to a Keplerian are uncorrelated noise, is invalid, and re- 
drawing them should worsen the quality of the Keplerian (s) plus trend fit, since scrambling 
will remove evidence of the trend. Thus, the fraction of these sets with Axl ^^^^ than that 
of the proper, unscrambled residuals, provides a measurement of the false alarm probability 
that the residuals to a Keplerian-only fit are correlated. 



3.2. 



Velocity Trends and Additional Companions in Known Exoplanet 

Systems 



The Catalog of Nearby Exoplanets (IButler et al.ll2006l ) contains 172 substellar compan- 
ions with msini < 24Mj^p orbiting 148 stars within 20 pc. Since then , at least 3 more 
systems have been announced, includin g a triple-Neptune (ILovis et al.ll2006f) . and two single- 
planet detections ( Johnson et al. 2006 ). ( Hatzes et al. 20061 ). Of these 151 systems, 24 show 
significant trends in addition to the Keplerian curve s of the known exoplanets. We have 
reanalyzed the radial velocities of iButler et al.l (120061 ) to determine the significance of these 
trends and to find evidence for additional trends using the FAP test des cribed in ^I3.1 [ Note 
that we have obtained additional data for some of these systems since Butler et al.l ( 2006 ) 
went to press. 



We confirm here 21 of the 24 trends in iButler et al.l (120061 ) to have FAPs below 1% (2 
others are in systems on which we have no data to test, and the third is HD 1196 4, discussed 



i n ^ 16.31) . We also confirm the trend in the 14 Her system, first announced in Naef et al. 



(120041 ) and analyzed more thoro ughlv in iGozdziewski. Konacki. fc Maciejewski J2006h. and 



^ 16.21 W e confirm the finding of lEndl et al.l (120061 ) that the trend reported in 



Marcy et al. 



(l2005bl ) for HD 45350 b is not significant (FAP=0.6 and xt increases with the introduction 
of a trend). 

We announce here the detection of statistically significant linear trends (FAP < 1%) 
around 4 stars already known to harbor a single exoplanet: HD 83443, GJ 436 (= HIP 57087), 
HD 102117, and HD 195019. GJ 436 will be discussed more thoroughly in an upcoming work. 
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(Maness et al., 2006, submitted) ). In one additional case, HD 168443, we detect a radial 
velocity trend with FAP < 1% in a system already known to have two exoplanets, indicating 
that a third, long-period companion may exist. We present the updated orbital solutions in 
Table E) 

HD 49674 has an FAP for an additional trend of ~ 2%, which is of borderline significance 
when we account for the size of our sample: we should expect that around 2 of our 100 
systems will prove to have FAFs ~ 2% purely by chance, and not because of an additional 
companion. We include the fit for HD 49674 with a trend in Table |31 but note here the 
weakness of the detection. 



4. Constraining Long Period Companions 

4.1. The Problem of Incomplete Orbits 

It is difficult to properly characterize the orbit of an exoplanet when the data do not 
span at least one complete revolution. After one witnesses a complete orbit of the planet 
in a single-planet system, subsequent orbits should have exactly the same shape (absent 
strong planet-planet interactions), and so one can interpret deviations as the effects of an 
additional companion. Before witnessing one complete orbit, one can easily misinterpret 
the signature of an additional companion as it is absorbed into the orbital solution for the 
primary companion. Even when only one planet is present, small portions of single Keplerian 
curves can easily mimic portions of other Keplerians with very different orbital elements. 



Table 2. Properties of Three Stars Hosting New Substellar Companions 



HD 


Hip # 


RA 


Dec. 


B-V 


V 


Distance 




log 9 


[Fe/H] 


V sin i 


Mass 


S 


AMy 


jitt 






(J2000) 


{J2000) 






(pc) 


(K) 


(cm s~ ^ ) 




(ms-1) 


(Mq) 




(m s 


24040 


17960 


03 50 22.968 


+17 28 34.92 


0.65 


7.50 


46.5(2.2) 


5853(44) 


4.361(70) 


0.206(30) 


2.39(50) 


1.18 


0.15 


0.65 


5. 


154345 


83389 


17 02 36.404 


+47 04 54.77 


0.73 


6.76 


18.06(18) 


5468(44) 


4.537(70) 


-0.105(30) 


1.21(50) 


0.88 


0.18 


-0.21 


5. 




14810 


03 11 14.230 


+21 05 50.49 


0.78 


8.52 


52.9(4.1) 


5485(44) 


4.300(70) 


0.231(30) 


0.50(50) 


0.99 


0.16 


0.64 


3. 



Note. — For succinctness, wc express uncertainties using parenthetical notation, where the least significant digit of the uncertainty, in parentheses, and that of the qu; 
are understood to have the same place value. Thus, "0. 100(20)" indicates "0.10 ± 0.020", "1.0(2.0)" indicates "1.0 ± 2.0", an d "1(20)" indicates "1 ib 2 0" . 
Data from columns 3—5, and 12 are from Hipparcos jPerryman fc ESA|[T99^^. columns 6—10 are from the SPOCS catalog jValenti fc Fischeil l2005bD . column 11 is 
IWright et al.| j2004| ) . and column 13 was derived from using the formula iiT IWrightl l |2005|) . Columns 6—10 for HIP 14810 were derived with the same methods used fc 
SPOCS catalog. 
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4.2. Constraining msinz and P 



When an RV curve shows significant curvature, it may be poss i ble to constrain the 
minimum mass (m sin i) and orbital period of the companion. iBrownl (120041 ) discussed the 
problem extensively, and IWittenmyer et al.l (120061 ) studied the significance of non-detections 
in the McDonald Observatory planet search without assuming circular orbits by injecting 
artificial RV signal s into program data to determine the strength of a just-recoverable signal. 
Wittenmyer et al.l (120061 ) reasonably assigned a broad range of eccentricities, < e < 0.6, 
with an upper limit they j ustified by the fact that over 90% of all known exoplanets have 
e < 0.6 JSutler et al.ll2006h . The presence of this upper limit greatly limits the number of 
pathological solutions to a given RV set. Below, we explore the nature of limits on mass and 
period implied by a given data set and how constraining e can improve those limits. 

Since msinz, not is the astrophysically interesting quantity in exoplanet detection, 
it is useful to transform into P, e, msini coordinates when considering constraints. The 
minimum mass of a companion can be calculated from the mass function, f{f^)i and the 
stellar mass, according to the relation: 



/(m) 



w? sin^ i 
(m + M,)2 



27rG' 



(1) 



where, in the minimum mass case (where sini = 1), we set m equal to msini. This relation 
allows us to fit for the minimum mass (which we refer to as m sin i for brevity) eliminating 
the orbital parameter K. 

Using Eq. [T], we can find the best-fit Keplerian RV curve across P — msini space, 
allowing e, u;, and 7 (the RV zero point) to vary at many fixed values of P and msini to 
map x^. 



Table 3. Updated Orbital Fits for 9 Exoplanets 



Planet 




Per 


K 


c 






trend 


m sin i 


a 


r.m.s. 




Nobs 






(d) 


(m s-1) 




(dog) 


(JD-2440000) 


(m/s/yr) 


(Mjup) 


(AU) 


(ms-1) 






HIP 14810 


b 


6.6742(20) 


428.3(3.0) 


0.1470(60) 


158.6(2.0) 


13694.588(40) 




3.91(55) 


0.0692(40) 


5.1 


1.4 


30 


HIP 14810 


c 


95.2914(20) 


37.4(3.0) 


0.4088(60) 


354.2(2.0) 


13679.575(40) 




0.76(12) 


0.407(23) 


3.3 


0.77 


21 


HD 49674 


b 


4.9437(23) 


13.7(2.1) 


0.29(15) 


283 


11882.90(86) 


2.6(1.1) 


0.115(16) 


0.0580(33) 


4.4 


0.62 


39 


HD 83443 


b 


2.985625(60) 


56.4(1.4) 


0.008(25) 


24 


11211.04(82) 


2.40(79) 


0.400(34) 


0.0406(23) 


8.2 


0.93 


51 


GJ 436 


b 


2.643859(74) 


18.35(80) 


0.145(52) 


353(24) 


11551.72(12) 


1.42(35) 


0.0682(63) 


0.0278(16) 


4.2 


0.93 


60 


HD 102117 b 


20.8079(55) 


11.91(77) 


0.106(70) 


283 


10942.9(3.0) 


-0.91(26) 


0.172(18) 


0.1532(88) 


3.3 


0.83 


45 


HD 168443 


b 


58.11289(86) 


475.9(1.6) 


0.5286(32) 


172.87(94) 


10047.387(34) 




8.02(65) 


0.300(17) 


4.1 


0.97 


109 


HD 168443 


c 


1749.5(2.4) 


298.0(1.2) 


0.2125(15) 


65.07(21) 


10273.0(4.6) 




18.1(1.5) 


2.91(17) 


4.1 


0.97 


109 


HD 195019 


b 


18.20163(40) 


272.3(1.4) 


0.0140(44) 


222(20) 


11015.0(1.2) 


1.31(51) 


3.70(30) 


0.1388(80) 


16 


1.5 


154 



Note. — For succinctness, we express uncertainties using parenthetical notation, where the least significant digit of the uncertainty, in parentheses, and that 
of the quantity are understood to have the same place value. Thus, "0.100(20)" indicates "0.100 zb 0.020", "1.0(2.0)" indicates "1.0 zb 2.0", and "1(20)" indicates 
"1 ± 20". 
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5. Two New Substellar Companions with Incomplete Orbits 

Here we consider HD 24040, a metal-rich ([Fe/H] = 0.21) GO V star at 46 pc (stellar 
characteristics summarized in Table [2]). This star shows RV variations consistent with a 
planetary companion with P ~ 15 y and msini ~ '''^Jup (Figure [2]), although longer 
orbital periods and minimum masses as high as msini ~ 30 cannot be ruled out. The RV 
data for HD 24040 appear in Table H 

Here and in §6.21 we use ^is a merit function and infer parameters for acceptable 
fits from increases of this function by 1, 4 and 9, which correspond to 1-, 2-, and 3-sigma 
confidence levels for systems with Gaussian noise. Because stellar jitter provides a source of 
pseudo-random noise which may vary on a stellar rotation timescale, the noise in RV residuals 
may be non-Gaussian. Thus, to the degree that the RV residuals are non-Gauassian, the 
translation of these confidence limits into precise probabilities is not straightforward. 
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Tabic 4. RV Data for HD 24040 



Time 


Radial Velocity 


Unc. 


(^JJJ-z44UUUU ) 


(m/ s) 


(m/ s) 


10838.773206 


-37.1 


1.4 


11043.119653 


-25.8 


1.5 


11072.039039 


-24.3 


1.4 


11073.002315 


-26.8 


1.2 


11170.876921 


-9.9 


1.4 


11411.092975 


3.6 


1.6 


11550.824005 


16.4 


1.3 


11551.863449 


16.5 


1.2 


11793.136725 


42.2 


1.2 


11899.945741 


50.7 


1.1 


12516.065405 


74.4 


1.3 


12575.951921 


87.6 


1.5 


12854.115278 


81.4 


1.2 


12856.115671 


82.1 


1.1 


13071.741674 


59.4 


1.2 


13072.799190 


55.6 


1.3 


13196.130000 


54.2 


1.3 


13207.120116 


57.3 


1.3 


13208.125625 


54.0 


1.2 


13241.091852 


57.5 


1.2 


13302.947025 


49.9 


1.3 


13339.945972 


46.1 


1.2 


13368.865139 


48.4 


1.2 


13426.809792 


36.5 


1.1 


13696.904468 


18.0 


1.1 


13982.061586 


6.8 


1.2 
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In this case, we have enough RV information to put an upper hmit on msini. Figure [3] 
shows foi' best-fit orbits in the P — msini plan. Fits with P as low as 10 y and msini 
as low as 5 Mj^p are allowed. Interestingly, the data (following the middle, = X^j^ + 4 
contour) exclude orbits with msini > SOMj^p, providing a "maximum minimum- mass" . 
Since without an assumption for the eccentricity (which we will make below), we cannot 
exclude orbits with msini as high as 30 Mj^p, there is a chance that this companion to HD 
24040 is a brown dwarf, or even stellar. 

A similar case is HD 154345, a G8 V star at 18pc (stellar characteristics summarized in 
Table [2]). This star shows RV variations remarkably similar to those of HD 24040 (Figure H]), 
but with an amplitude about 6 times smaller. In this case, the maximum msini is only 
around 10 Mj^p, giving us confidence that this object is likely a true exoplanet, and masses 
as low as 1 Mj^p are allowed. The RV data for HD 154345 are in Table [5l We summarize 
the orbital constraints for these objects in Table [61 
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Fig. 2. — RV curve for HD 24040 with data from Keck. The best- fit Keplerian is poorly 
constrained due to incomplete coverage of the orbit. The fit shown here is for P — 16.5 y 
and msini = 6.9 Mt„„, one of a family of adequate solutions. 
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Fig. 3. — Contours of and e in P — msini space of best-fit orbits to the RV data of HD 
24040 (Fig. [2]), with in grayscale. The sohd contours mark the levels where increases 
by 1, 4 and 9 from the minimum. The dashed contours mark levels of the eccentricity of 
0.2, 0.6, and 0.9, Planets with e > 0.6 are rare, implying that this object is unlikely to have 
a period longer than 100 y. The orbit is largely unconstrained, but msini has a maximum 
value around 20 Mj^p for orbits with e < 0.6. The position of the cross at 16.5 y and 6.9 
M T„„ represents the solution plotted in Fig. [2], one in a family of adequate orbital solutions. 
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Table 5. RV Data for HD 154345 



Time Radial Velocity Unc. 
(JD-2440000) (m/s) (m/s) 



10547.110035 


6.8 


1.4 


10603.955845 


8.6 


1.4 


10956.015625 


15.1 


1.5 


10982.963634 


13.1 


1.4 


11013.868657 


16.3 


1.5 


11311.065486 


16.6 


1.6 


11368.789491 


18.7 


1.5 


11441.713877 


23.0 


1.4 


11705.917836 


30.3 


1.5 


12003.078183 


30.4 


2.3 


12098.916539 


37.1 


1.5 


12128.797813 


34.7 


1.7 


12333.173299 


38.1 


1.6 


12487.860197 


35.5 


1.6 


12776.985463 


29.9 


1.6 


12806.951852 


19.5 


1.6 


12833.801030 


27.3 


1.4 


12848.772037 


25.3 


1.5 


12897.776562 


26.5 


1.5 


13072.046921 


18.9 


1.6 


13074.077766 


21.2 


1.4 


13077.128090 


20.5 


1.4 


13153.943171 


15.2 


1.6 


13179.992454 


20.6 


1.5 


13195.819190 


17.5 


1.4 


13428.162502 


10.09 


0.78 


13547.914433 


11.44 


0.80 


13604.829999 


6.08 


0.78 


13777.155347 


7.5 


1.5 


13807.077257 


2.4 


1.4 
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We can put more stringent co nstraints on these orbits by noting that ninety percent of 
all known exoplanets have e < 0.6 JSutler et allbooeh . For both HD 24040 b and HD 154345 
b, the high-peri od solutions all have hig h eccentricities (the dashed contours in Figs. [3] and 
EJ. If, following I Wittenmyer et al.l (120061 ) . we therefore assume that e < 0.6 for these objects, 



and that the true values of m sin i and P lie within the limits of the middle (x^ 



mm 



+ 4) 



contour, then we can constrain 5Mj^p < msini < 20Mj^^p and 10 < P < 100 y for HD 
24040 b, and 0.8 Mj^p < msini < 10 Mj^p and 7 < P < 100 y for HD 154345 b. 



6. Mining Velocity Residuals for Additional Exoplanets 

6.1. Velocity Residuals Suggesting Additional Companions 

For exoplanetary systems in which an additional, low-amplitude signal is not well- 
characterized by just a linear trend — for instance, where there is significant curvature 
(e.g. HD 13445 or HD 68988) or even multiple orbits (e.g. GJ 876) — a full, multi-planet fit 
is needed to properly characterize the system. In this case, we can apply an FAP analysis 
similar to the one in § 13.11 testing the (A^ + l)-planet hypothesis versus the null hypothesis 
of planets plus a trend plus noise, where N is the number of previously confirmed plan- 
ets. This is a much more computationally intensive procedure than that of § 13.11 since we 
are introducing 5 new, non-linear, highly covariant parameters (P, e,a;,Tp, and K), so we 
have performed only 50-100 trials. In most cases the low-amplitude signal we seek is much 
weaker than that of the known planet (s). This means we have good initial guesses for the or- 
bital parameters of the established exoplanets, and that those parameters are usually rather 
insensitive to those of the additional companion, easing the difficulty of the simultaneous 
11-parameter fit (16-parameter for existing double systems). 

As in § 13.11 we calculate the improvement in the goodness-of-fit parameter, Axl = 

1 ■ — 1 , , , i with the introduction of an additional exoplanetary com- 
^u,N+i planets ^!^,Afplanets+trend ^ 

Table 5 — Continued 



Time Radial Velocity Unc. 
(JD-2440000) (m/s) (m/s) 



13931.955714 
13932.913019 



1.33 0.72 
1.98 0.70 
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Fig. 4. — RV curve for HD 154345 with data from Keck. The best-fit Keplerian is poorly 
constrained due to incomplete coverage of the orbit. The fit shown here is for P = 35.8 y 
and msini = 2.2 Mj„„ (marked in Fig. [5]), one in a family of adequate orbital solutions. 
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Fig. 5. — Contours of and e in P — msini space of best-fit orbits to the RV data of HD 
154345 (Fig. H]), with in grayscale. The sohd contours mark the levels where increases 
by 1, 4 and 9 from the minimum. The dashed contours mark levels of the eccentricity of 
0.2, 0.6, and 0.9. Planets with e > 0.6 are rare, implying that this exoplanet is unlikely 
to have a period longer than 100 y. The orbit is largely unconstrained, but msini has a 
maximum value around 10 Mj^p for orbits with e < 0.6. The white cross at 36 y and 2.2 
Mj^p represents the solution shown in Fig. HI 
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panion compared to a fit including only an additional trend. We compare this reduction to 
that of mock data sets bootstrapped as the sum of the best-fit solution with a trend, plus 
noise, drawn, with replacement, from the residuals of the actual data to this fit. We then 
construct an FAP as the fraction of mock sets that saw a greater reduction in xl with the 
introduction of an additional planet than the genuine data set. 

A low FAP for the presence of a second planet is not tantamount to the detection of 
an additional exoplanet. It is only a sign that the null hypothesis is unlikely, i.e. that the 
distribution of residuals is not representative of the actual noise in the system or that the 
presumed orbital solution from which the residuals were drawn is in error. This would be 
the case if, for instance, if the residuals are correlated due to non-Keplerian RV variations 
(such as systematic errors or astrophysical jitter). 

The fits discussed here are purely Keplerian and not dynamical. In particular, fits which 
produce unstable or unphysical orbits are allowed. More sophisticated, Newtonian fits (e.g. 



Rivera et al.ll2005l ) would better constrain the orbits of multiple planet systems. 



6.2. Long-Period Companions with Incomplete Orbits 



We have identified 8 other systems in which the FAP for an additional Keplerian vs. a 
simple trend is below 2%: HD 142, HD 13445, HD 68988, 23 Lib (= HD 134987), 14 Her, 
T Boo (= HD 120136), HD 183263, and HD 187123. In addition, we have identified a ninth 
system, HD 114783, which has a compelling second Keplerian despite a slightly larger FAP 
(6%). We summarize the orbital constraints for these objects in Table O 

— HD 142: Most of the RV data for HD 142 show a sim ple linear trend superimposed 
on the known K = 34 m/s, 350-d orbit (ITinney et al.ll2002l ). HD 142 is known to have a 
stellar companion {V = 10) ( Poveda et all 1994 ). which could explain the trend. The first 
two data points, taken in 1998-9, are significantly low, producing a low FAP for curvature 
(< 1%). HD 142 has B — V= 0.52, indicating it is a late F or early G star, suggesting it 



may have moderate jitter (~ 5 ms ^, IWrightl (j2005l )). so we therefore view the low FAP 



for curvature, apparently based on only two low points, with suspicion. If the curvature is 
real, it is consistent with an exoplanet with period longer than the span of the observations 
(P > 10 y) with a minimum mass of at least 4 Mj^p. 



HD 13445 (= GL 86) has a known planet with P = 15.76 d flQueloz et al.l 12000 



Butler et all l2006h .Superimposed on that Keplerian velocity curve is a velocity trend of 
roughly -94 m/s/yr during the past 9 yea r s, apparently consistent with t he brown dwarf cora - 

Els et al.l boOW . Ichauvin et al.l J2006h : lOueloz et al.l J200oh . 



panion previously reported by 
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1998 2000 2002 2004 2006 

Date 

Fig. 6. — RV curve for HD 142, a multiple companion system, with data from from AAT, 
. The previously known inner planet has P = 350 d, and the outer companion is poorly 
constrained, but consistent with the known stellar companion. The data are inconsistent 
with a linear trend, mostly because of the first two data points. 
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stellar or brown-dwarf companion (iQueloz et al.ll2000l ). There is a a hint of curvature in 
these residuals to the inner planet, but not enough to put meaningful constraints on this 
outer object beyond that fact that its period is longer than the span of the observations 
(~ 10 y) and msmi> 22 Mj^p. 

— HD 68988 shows definite signs of curvature in the residuals to the 1.8 Mj^p inner 
planet (as Fig. [7] shows). Fig. [8] shows the outer companion has msini < 30Mj^p, and the 
assumption of e < 0.6, using the middle contour, further restricts msini < 20Mj^p, and 
P < 60y. 

— HD 114783 shows curvature in its residuals, and may have experienced both an RV 
minimum (in 2000) and maximum (in 2006) as the RV curve in Fig. [9] shows. The data are 
only moderately inconsistent with a linear trend, however (FAP = 6%), indicating that the 
outer companion's orbit is still underconstrained. 

— 23 Lib (= HD 134987) shows signs of curvature in the residuals to the known inner 
planet . The signal appears as a change in the level of otherwise fiat residuals between 2000 
and 2002 of 15 m/s (see Fig. (TUj). This suggests an outer planet on a rather eccentric orbit 
which reached periastron in 2001. The small magnitude of this change in RV suggests a 
low-mass object, but the incomplete nature of this orbit makes us less than certain that it 
is due to an exoplanet. 



— 14 H er (= HD 145645): This star has a known tren d (INaef et al.ll2004l ) and has been 
analyzed by iGozdziewski. Konacki. fc Maciejewskil (120061 ) as a possible resonant multiple 
system. The previously-known planet has msini = 4.9 Mj^p and P = 4.8 y, but the 
character of the second companion is uncertai n. Combining o ur data with the published 
ELODIE data from the Geneva Planet Search (INaef et al.l 120041 ) (Fig. [TT]) provides a good 
picture of the system. The character of the orbit of the outer planet is unconstrained, and 
several equally acceptable but qualitatively distinct solutions exist. One is a long period, 
nearly circular orbit like the one shown in Fig. 10 and a mass near 2 Mj^p. Other solutions 
include a 3:1 resonance with the inner planet. The next few years of observation should 
break this degeneracy. The degeneracy may also be broken by high contrast, high resolution 
imaging, and we suggest that such attempts be made on this interesting system. 

— HD 183263 shows definite signs of curvature in the residuals to the known inner 
planet (as Fig. [12] shows), but too little to constrain the mass of the distant companion. 
Fig. [12] shows that there is little meaningful constraint on the orbit beyond P > 7 j and 
msini > 4Mj^p. Even the assumption of e < 0.6 allows for msini > 13, so the planetary 
nature of the companion is very uncertain. 



r Boo (= HD 120136) has residuals to the fit for the known inner planet which show 
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Fig. 7. — RV curve for HD 68988, a multiple-companion system, with data from Keck. 
The previously known inner planet has P — 6.28 d, and the outer companion is poorly 
constrained, but likely has msini < 20Mj^p and P < 60y. 
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Fig. 8. — Contours of and Cc in — (msini)c space for the best double- Keplerian fits to 
the RV data of HD 68988 (Fig. [7j), with in grayscale. The solid contours mark the levels 
where increases by 1, 4 and 9 from the minimum. The dashed contours mark levels of 
the eccentricity of 0.2, 0.6, and 0.9. Assuming e < 0.6, we can constrain 6Mj^p < msini < 



20 M 



Jup 



and 11 < P < 60 y. 
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Fig. 9. — RV curve for HD 114783, a multiple-companion system, with data from Keck. 
The previously known inner planet has P = 495 d and msini = l.lMj^p. The residuals 
are only moderately inconsistent with a linear trend (FAP = 6%), indicating that the outer 
companion is poorly constrained. 
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Fig. 10. — RV curve for 23 Lib (= HD 134987), a multiple-companion system, with data 
from Keck and AAT. The previously known inner planet has P = 258 d and msini = 1.62. 
The outer companion is poorly constrained. The orbital parameters of the inner planet are 
not significantly changed with a two-parameter fit. 
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1994 1996 1998 2000 2002 2004 2006 




Fig. 11. — RV curve for 14 Her (= HD 145675), a system with multiple companions. Crosses 
repre sent data from t he ELODIE instrument operated by the Geneva Planet Search (taken 
from lNaef et al.ll2004l ). and large filled circl es represent data t aken at Keck Observatory by 
the California and Carnegie Planet Search (IButler et al.ll2006l ). Error bars represent quoted 
errors on individual velocities; for some points the error bars are smaller than the plotted 
points. The combined data set shows a long-period companion with P > 12y and msini > 5 



M 



jyp. The previously known inner planet has P = 4.8 y and msini = 4.9 Mj^p. 
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Fig. 12. — RV curve for HD 183263, a multiple-companion system, with data from Keck. 
The previously known inner planet has P — 635 d and msini = 3.8Mj^p, and the outer 
companion is poorly constrained. 
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Period (yr) 

Fig. 13. — Contours of and Cc in Pc — (msini)c space for the best double-Keplerian fits to 
the RV data of HD 183263 (Fig. [T2l) . with in grayscale. The sohd contours mark the levels 
where increases by 1, 4 and 9 from the minimum. The dashed contours mark levels of the 
eccentricity of 0.2, 0.6, and 0.9. P and msini for this companion are poorly constrained. 
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evide nce of a long-period companion which have been discussed elsewhere (IButler et al. 



20061 ). Analysis of the distant companion is complicated by the lower quality of the data 



during the apparent periastron in 1990. The current best fit suggests a period greater than 
15 years, but is otherwise unconstrained. iPoveda et al.l (119941 ) report that r Boo has a faint 
(V=10.3) companion (sep. 5.4") which may be the source of the RV residuals. 



HD 187123 is known to host a 0.5 Mj^p "Hot Jupiter" in a 3 day orbit (IButler et al. 



19981 ). Observations over the subsequent 8 year s have revealed a trend of -7.3 m/s in the 



residuals to a one planet fit (IButler et al.ll2006l ). In 2001, the trend began to show signs 
of curvature, and in 2006 a it became clear that the residuals had passed through an RV 
minimum (see Fig. [T^ . Fig. [15] shows the ^^^id e contours in P — msini space. In 
this case, the e = 0.6 contour and middle contour provide the following constraints: 
2 Mj^p < m sin i < 5 Mj^p and 10 < P < 40 y. 



6.3. Short-Period Companions to Stars with Knovi^n Planets 

We now consider known single-planet systems with low FAPs for second planets whose 
best-fit solutions have periods shorter than the span of observations. We have identified 
eight such systems, and we discuss them below. 

Five stars appear to exhibit coherent residuals (FAP < 2%) to a one planet fit, but in 
all cases the best two-Keplerian fits are not compelling (as noted in § 13. 2[ in a sample of 
100 known planet-bearing stars, we expect around 2 to exhibit residuals coherent at this 
level purely by chance). These possible companions do not appear in Table [3] because the 
tentative nature of these signals do not warrant publication of a full orbital solution with 
errors. 



HD 11964 was announced in iButler et al.l (120061 ) as having a planet with a 5.5 y 



orbital period and a linear trend. We find an FAP for a linear trend to be 6%, suggesting 
that while the inner planet is real, the trend is not. We find an FAP for a second planet 
to be < 2%, and a best-fit solution finds an inner planet with P = 37. 9d. This star sits 2 
magnitudes above the main sequence, and the residuals to t he known planet are consistent 



with the typical jitter for subgiants of 5.7 m/s (jWrightll2005l ). so this signal could represent 
some sort of correlated noise. This very low amplitude signal {K = 5.6 m/s) will thus require 
much more data for confirmation. 

— HD 177830 is already known to have a Jupiter- mass object in a nearly circular, 1.12 
y orbit. This remarkable system has a low FAP < 1% for a second planet vs. a trend. Two 
good two-planet solutions exist for this system: the first has P = 111 d and msinz = 0.19 
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Fig. 14.— RV curve for HD 187123, with data from Keck, showing the 0.5 Mj^p"Hot 
Jupiter" and the outer companion of uncertain period and mass. 
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Fig. 15. — Contours of and Cc in Pc — sin i)c space for the best two-planet fits to the 
RV data of HD 187123 (Fig. [HI), with iii grayscale. The solid contours mark the levels 
where increases by 1, 4 and 9 from the minimum. The dashed contours mark levels of 
the eccentricity of 0.2, 0.6, and 0.9. Planets with e > 0.6 are rare, implying that this object 
is unlikely to have a period longer than 40 y or msinz greater than 5 Mt,,„. 
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j^p, the second has P = 46.8 and msinz = 0.16 Mj^p. 
magnitudes above the main sequence, and the residuals to t he known planet are consistent 
with the typical jitter for subgiants of 5.7 m/s (IWrightll2005l ). so this signal could represent 
some sort of correlated noise. 

— 70 Vir (= HD 117176) is a subgiant with a massive, 116.6 d planet on an eccentric 
orbit (e = 0.39). The FAP for a second planet is 2%, but the best-fit second planet is not 
persuasive: P = 9.58 d, and K = 7 m/s. The typical internal errors for this target are 5.4 
m/s, making a bona fide detection of a 7 m/s planet very difficult. We suspect that this 
signal is an artifact of stellar jitter, possibly due to the advanced evolution of the star. 

— HD 164922 has a known planet with a 3.1 y orbital period. For this star, the FAP 
for a second planet is < 1%. The best fit for this second planet has P = 75.8 d and 
msinz = 0.06Mj^p. The amplitude of this signal is extremely low — only K = 3 m/s — 
making this an intriguing but marginal detection. 

— HD 210277 is already known to host a planet with a 1.2 y orbit. The FAP for a 
second planet is 2%, and the best-fit second Keplerian has K = 3 m/s signal and P = 3.14 
d, and a 2% FAP. The best-fit orbit has e = 0.5, which is unlikely given that nearly all known 
Hot Jupiters have e < 0.1 (although the presence of the 1.2 y, e = 0.5 outer planet could 
be responsible, in principle, for pumping an inner planet's eccentricity.) The extremely low 
amplitude of this planet makes the exoplanetary nature of this signal very uncertain. 

T hree additional stars w ith low FAPs are of a very e arly spectral type (F7-8): HD 



This star sits more than 3.5 



89744 flKorzennik et al.ll2000f ). HD 108147 (IPepe et al.ll2002f ). and HD 208 487 flTinnev et al. 



20051 ) . Their low activity yields a low jitter in the estimation of IWrightl (120051 ). but this is 
likely underestimated due to poor statistics: the California and Carnegie Planet Search has 
very few stars of this spectral type from which to estimate the jitter. For HD 89744 and HD 
108147 we suspect that, the low FAP of < 2% is an artifact of coherent noise, since in our 
judgment neither case shows a compelling evidence of a second Keplerian of any period. 

HD 208487 has a very low FAP (< 1%) despite the modest r.m.s of the residuals to a one 
planet fit of 8 m/s. We suspect that s tellar jitter is the likely source of these variations. This 
star was discussed by iGregoryl (120051 ) . who applied a Bayesian analysis to the published RV 

Jup' 



data, concluding that a second plan et was likely, having P = 998^^^ d and 



0.5M 



Gozdziewski fc Migaszewskil (120061 ) also studied the published data, and suggested a planet 
with P = 14.5 d. We note here two plausible solutions apparent i n our data. The first, with 
P ~ 1000 d and msini ~ 0.5, is consistent with the solution of iGregoryl (120051 ). We also 
fi nd, however, an additional so l ution of equal quality with P = 28.6 d (double the period 



of Gozdziewski fc Migaszewski ( 20061 )) and msini = 0.14Mj^p. This second solution has a 
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period uncomfortably close to that of the lunar cycle (we often see this period in the window 
function of our observations due to our tendency to observe during bright time). For both 
solutions K = 10 m/s. We reiterate that we feel that the early spectral type of this star 
alone can account for the observed RV residuals. 



Gozdziewski fc Migaszewskil (120061 ) analyzed our published RV data and found a low 



FAP for the existence of a second planet in orbit around HD 188015. Their FAP, however, 
is measured against a null hypothesis of a single Keplerian plus noise, thus ignoring the 
linear trend. We find that Axl, the improvement of the goodness-of-fit parameter with the 
introduction of a second Keplerian versus a trend to be very small — in fact 60% of our 
mock data sets showed greater improvement. We therefore find no motivation to hypothesize 
th e existence o f an add itional, short-period planet; the single planet and trend announced 
in 



Marcy et al.l (l2005bl ) are sufficient to explain the data. 



Gozdziewski fc Migaszewskil (120061 ) also found a low FAP for a se cond planet i n orbi t 



about HD 114729, with a period of 13.8 d. Using the data set from iButler et al.l (120061) 



which contains 3 recent RV measurements taken since the publication of iButler et al. 



(their source of RV data), we find no such signal, and a large FAP for a second planet. We 
suspect our results may differ because the additional data provide for a slightly better fit to 
the known exoplanet, changing the character of the residuals and destroying the coherence 
of the spurious 13.8 d signal. 



7. HD 150706 

HD 150706 b, a purported l.OMj^p eccentric planet at 0.8 AU, was announced by the 
Geneva Extrasolar Planet Search Team (2002, Washingtor i conf erence "Scientifi c Frontiers in 



Research in Extrasolar Planets"; Udry, Mayor & Queloz. |2003| ) and appears in iButler et al. 



( 120061 ): however, there is no refereed discovery paper giving details. 



We have made eight precision velocity measurements at Keck observatory from 2002 
through 2006. These velocities show an RMS scatter of 12.1 m/s, inconsistent with the 
reported 33m/s semi-amplitude of HD 150706 b. The RMS to a linear fit is 8 m/s, which is 
adequately explained by the expected jitter for a young (700±300 Myr) and active star like 
HD 150706. We therefore doubt the existence of a 1.0 Mj^peccentric planet orbiting HD 
150706 at 0.8 AU. 
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8. Discussion 



As noted in § I3.2[ prior to tliis work 24 of tlie 150 nearby stars known to liost exoplanets 
(including 14 Her and excluding HD 150706) show significant trends in their residuals and 
19 host well-characterized multiple planet systems. One of these trends is likely spurious 
(HD 11964), and at least 3 others may be due to stellar or brown dwarf companions (HD 
142, HD 13445, and r Boo). We have announced here the detection of an additional 5 
trends for known planet-bearing stars, 2 new single sys tems, and on e new multiple system 



(HIP 14810, which appears as a single-planet system in iButler et al.l (120061 )). We have also 
confirmed that the previously announced trends for HD 68988 and HD 187123 are likely due 
to planetary-mass objects. This brings the total number of stars with RV trends possibly 
due to planets to 22, the number of known multiple-planet systems to 22, and the number 
of nearby planet-bearing stars to 152. This means that 30% of known exoplanet systems 
show significant evidence of multiplicity. Considering that t he mass distribution of planets 



increases steeply toward lower masses (IMarcy et al.l l2005al ) , our incompleteness must be 



considerable between 1.0 and 0.1 Jupiter- masses. Thus, the actual occurrence of multiple 
planets among stars having one known planet must be considerably greater than 30%. 

From an anthropocentric perspective, this frequency of multiplicity suggests that in 
some respects, the Solar System is not such an aberration. Our Sun has 4 giant planets, and 
it appears that such multiplicity is not uncommon, although circular orbits are. 

From a planet-hunting perspective this result is quite welcome as well, since it means 
that the immediate future of RV planet searches looks bright. As our temporal baseline 
expands, we will become sensitive to longer-period planets. Our search is just becoming 
sensitive to true Jupiter analogs with 12 year orbits and 12 m/s amplitudes. A true Saturn 
analog would require 15 more years of observation. As our precision improves we will become 
sensitive to lower-mass planets, which may be the richest domain for planets yet. 
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Table 6. Mass constraints for some substellar companions with incomplete orbits 



Object 




Per 


msini 


a 






(y) 


(Mjup) 


(AU) 


HD 24040 b 


10 


— 100 


5 — 20 


5 — 23 


HD 68988 c 


11 


— 60 


11 — 20 


5 — 7 


HD 154345 b 


7- 


- 100 


0.8 — 10 


4 — 25 


HD 187123 c 


10 


— 40 


2 — 5 


5 — 12 



Note. — These contraints correspond to the 
extrema of the are given by y'^- + 4 contour in 
P — msini space for orbits with e < 0.6. 



